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Abstract 

The AC loss in ReBCO coated conductors is large in situations when 
the conductors are subjected to a considerable magnetic field, like in ro- 
tating machines, transformers and high-field magnets. Roebel cables can 
reduce the AC loss in these cases. However, computer simulations are 
needed to interpret the experiments, understand the loss mechanisms, re- 
duce the AC loss by optimising the Roebel cable and design the cryogenic 
system. In this article, we simulate and discuss the AC loss due to an ap- 
plied magnetic field making an arbitrary angle with the cable and taking 
into account a realistic anisotropic field dependence of the critical cur- 
rent density. We study the AC loss in the superconductor parts for the 
limits of very high coupling currents and completely uncoupled strands. 
The simulations for the uncoupled case also describe a double pancake 
coil with no transport current. For the simulations, we use two different 
numerical methods with complementary strengths. This serves as a mu- 
tual check of the correctness of the simulation results, which agree to each 
other. Opposite than expected, we found that the AC loss does not only 
depend on the perpendicular component of the applied magnetic field. We 
also found that the AC loss for applied fields with an orientation below 7 
degrees with the strands surface is reduced more than one order of mag- 
nitude comparing to an untransposed cable. Therefore, we recommend to 
use Roebel cables for windings with important parallel components of the 
magnetic field, such as transformers and high-field magnets. 

1 Introduction 

ReBCO coated conductors present several advantages compared to other super- 
conductors, such as their higher operation temperature (up to the liquid nitrogen 
temperature) and their good performance in large magnetic fields. Nowadays, 
coated conductors are produced in long lengths with uniform properties, so they 
are ready to be used in market applications. In particular, power applications 
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are very promising, such as fault-current limiters, cables, rotating machines (in- 
cluding wind generators) and transformers. However, the AC loss in rotating 
machines and transformers is too high because of two reasons. First, the compo- 
nent of the AC magnetic field perpendicular to the tape is large. Second, these 
applications usually require high-current conductors made of several coated con- 
ductors. This results in thick cables with a substantial AC loss contribution 
from the parallel AC magnetic field. A similar situation arises in high-field [T] 
or large-scale DC magnets (such as for particle accelerators P] or fusion). In 
this case, the it is still desirable to reduce the AC loss because it limits the 
ramp rate and generates significant histeresis loss, especially for liquid-helium 
cooled magnets. Roebel cables [3l El [5j [6j [7] reduce the ac loss caused by both 
components of an external magnetic field [H [5J |H] . In spite of the progress in 
Roebel cables, many aspects of their AC loss remain unknown [2]. A similar 
situation exists for pancake coils, widely used in windings. 

Computer simulations are needed to interpret the experiments, understand 
the loss mechanisms, reduce the AC loss by optimising the Roebel cable or 
pancake coil and design the cryogenic system. However, to the best of our 
knowledge, there are no published simulations for the angular dependence of 
the magnetisation AC loss in Roebel cables and pancake coils. 

The published simulations for the angular dependence of the magnetisation 
AC loss are only for a single tape or stacks of tapes. For a single tape, the pub- 
lished results are for a constant critical current density [5], Jc, or an anisotropic 
magnetic field-dependent Jc [111 [HI [HI [131 E] • However, the anisotropies and 
fielcl]J dependences utilized in these articles are not realistic, since they take into 
account only a generic dependence [TT], a dependence extracted from data at 
high applied fields (>0.5 T) [T0l[T5], or without the correction for the self-field 
[12l[T3]. Moreover, |9l [11] neither calculate the complete AC cycle nor the AC 
loss. The work in [T3] was for only two angles (90° and 15°). For stacks, there is 
only one work for the coupled case (allowing the magnetisation currents to close 
in different tapes) and relatively high applied field amplitudes (50 mT) [T3]. In 
addition, all these simulations were done assuming either a simplified anisotropy 
of Jc by an elliptical dependence or assuming an isotropic superconductor. In 
contrast, actual coated conductors present complex anisotropies [TOl [TTl fT51 [T^ . 

Simulations on Roebel cables are only for perpendicular applied fields [^[^ 
or transport currents [5D1 HH [52] . For double (or single) pancake coils (or stacks 
of tapes) the magnetisation AC loss of stacks of tapes was simulated only for 
perpendicular applied fields in [23 l [24 l [20 l [25 l [2T 1 |26] . 

The only measurements on the angular dependence of the AC loss are in [27| . 
However, there exists extensive work on coated conductors [551 [2S| and stacks 
of them jSni ISIl [32 ■ AH these measurements for single tapes and stacks were for 
coated conductors at an early stage of development, with poor artificial pinning 
and anisotropies approximately elliptical. As a consequence, they may not be 
representative for the present nano-engineered material [16l [171 [HI ■ For Roebel 
cables there are also measurements for perpendicular applied fields |^[351 [71[?7 ] . 
parallel ones [31] and transport currents HOI [33 . 

In this article, we simulate and discuss the AC loss due to an applied mag- 
netic field at an arbitrary angle with the cable and we take into account a 

^In order to avoid unnecessary repetition, we refer to the magnetic field as simply "field" 
except when the electric and magnetic fields could be confused with each other. 
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realistic anisotropic field dependence of the critical current density. We study 
the AC loss in the superconductor parts for the limits of very high coupling cur- 
rents (coupled case) and completely uncoupled strands (uncoupled case). The 
2D simulations for the uncoupled case also describe a double pancake coil with 
no transport current. Therefore, the results and discussions for the uncoupled 
case are also valid for double pancake coils. Additionally, we also discuss the 
details of the AC loss in a single tape. For the simulations, we use two differ- 
ent numerical methods with complementary strengths: the Minimum Magnetic 
Energy Variation (MMEV) and a Finite Element Method (FEM) with the H- 
formulation and edge elements (see section [O]) . This serves as a mutual check 
of the correctness of the simulation results. Moreover, it is also a check of the 
applicability of the sharp E{J) relation from the critical-state model {E and 
J are the electric field and the current density, respectively), because MMEV 
assumes this sharp E{J) relation and FEM, a smooth one. 

This article is structured as follows. In section [2] we outline the simulation 
models and the anisotropic field dependence of Jc for the calculations. In section 
13] we present and discuss the results for a single tape and a Roebel cable in the 
coupled and uncoupled cases. For the Roebel cable, we do not only present the 
AC loss but also the field and current distribution for some cases. Finally, in|3] 
we present our conclusions. 

2 Models 

In this section, we first outline both simulation methods, their complementary 
strengths and some technical details (section l2.ip . Afterwards, we detail the 
anisotropic field dependences of Jc in the simulations: two anisotropic (a 
realistic one for YBCO and an elliptical one) and one isotropic field dependences 
(section [221) ■ 

2.1 Simulation methods 

In this article, we use two different numerical methods to obtain the current dis- 
tribution, the magnetic field and the AC loss, as we did in a previous work j21| . 
These models are the Minimum Magnetic Energy Variation (MMEV) method 
and the Finite Element Method (FEM) with iJ-formulation and edge elements. 
The comparison of their results serves as a mutual check of the correctness of the 
methods, as well as of the assumption of the critical-state model for modelling 
high-temperature superconductors. 

These methods present different strengths. The MMEV method is generally 
faster than FEM [21]. This is an advantage to calculate situations with many 
tapes, such as windings [3S1[37]- Since MMEV is a user-programmed software 
(the program for this article is written in FORTRAN language), it is possible 
to control the processes, make further improvements and include it in larger 
programs. In addition, it is possible to exploit the vast collection of existing 
free-source numerical routines in FORTRAN or C-I--I-. The FEM model is more 
versatile because it uses a commercial software [38]. In particular, it can simu- 
late non-linear magnetic materials interacting with the superconductor [39] and 
multi-physics problems, as for example the coupling of electromagnetic and ther- 
mal effects [10] • Moreover, in contrast to MMEV, the FEM model can simulate 
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over-critical currents because it takes into account a smooth E{J) relation. 

Both methods are 2D models, reducing the problem to solving the cross- 
section of the cable. In more detail, they assume that the transposition length 
of the cable is much larger than its thickness and width, so the cable is well 
approximated by a set of infinitely long tapes parallel to each other, as we 
detailed in pT| . 

2.1.1 The Minimum Magnetic Energy Variation (MMEV) method 

The MMEV method assumes the sharp E{J) relation of the critical-state model, 
where E is the electrical field. It is based on a variational principle proposed by 
Prigozhin [41] , which finds the current distribution by minimising the magnetic 
energy variation, and a fast non-standard minimisation routine. This routine 
has been developed incrementally in several articles. First, Sanchez and Navau 
solved a cylinder in an applied magnetic field under certain restrictions [42) . 
Later on, Pardo et al. developed the general method for tapes under any com- 
bination of applied magnetic field and transport current [43| . The latest stage of 
the method is published in [44l |45l [ST] , where [44l |45] and [45] take into account 
the field dependence of Jc and the interaction with linear magnetic materials, 
respectively. Recently, we applied the method to Roebel cables [101 HI], but 
using a constant Jc- 

2.1.2 The Finite Element Method (FEM) model 

The FEM model assumes a smooth E{J) relation, usually a power law E{J) = 
Ec{J/ Jc)^ , where Ec is the voltage-per-length criterion for the critical current 
density Jc (usually set equal to 10~^ Vm~^) and n is the flux-creep exponent. 
In this work, we assume that n is constant, neglecting the magnetic field depen- 
dence of this parameter. The state variable of the FEM model are the magnetic 
field components; as a consequence the implementation of the Jc(B) dependence 
is straightforward since the magnetic flux components are immediately available 
from the state variables by means of the B = hqH relation. More details about 
the model implementation can be found in |46) . 

2.1.3 Tape and Roebel cable parameters 

The anisotropic field dependences of Jc for the simulations are detailed in 12.21 
These dependences are based on measurements [H] of YBCO coated conductors 
from Superpower, Inc. [47j . The functions for Jc describe the experimental 
anisotropic field dependence with several degrees of approximation: realistic, 
elliptic and isotropic. For all cases, we do not take into account the metal parts 
of the tapes in the calculations, ignoring the effects of eddy currents. 

For the single tape, the dimensions for the simulations are 4.16 mm for the 
width and 1.4 fiui for the thickness. 

For the Roebel cable, we chose the geometry of a cable composed of 14 
strands, which was manufactured at the Karlsruhe Institute of Technology ^U\ . 
The strands are 1.98 mm wide, their lateral separation is 200 /im and their ver- 
tical separation (i.e. the distance between the superconducting films) is 140 /im. 
The total critical current of the cable at 77 K (determined with the 1 /iV/cm 
criterion over a distance of 30 cm) is 465 A. In the simulations, the strands 
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(from now on we call them "tapes") have the following dimensions: width 
w = 1.98 mm and thickness d = 1.4 ^m. For the power-law resistivity used 
in the FEM simulations, Ec = 10^^ V/m and n = 35 (except stated otherwise). 
The frequency of the applied field is 100 Hz. 

2.1.4 Coupled and uncoupled cases 

For the cable simulations, we distinguish between the coupled and uncoupled 
cases. The coupled case assumes that the resistance per transposition length 
between the tapes is very small, so the magnetisation current loops can close 
freely between any tape of the cable. The uncoupled case assumes a very large 
resistance per transposition length between the tapes, so the current loops must 
close within each tape. From the computation point of view, this means that 
for the coupled case there is only one current constrain (zero net current in the 
whole cable), while for the uncoupled one there are as many current constrains 
as tapes (zero net current in each tape). 

A cable made of untransposcd tapes corresponds to the coupled case because 
the tapes are interconnected together at the current leads. Indeed, Polak et 
al. |48| experimentally found that for tape lengths of around 10 cm or above, 
interconnected tapes in parallel are already fully coupled. 

A double pancake coil exactly corresponds to the uncoupled case. Then, all 
the results and discussions for the Roebel cable in the uncoupled case are also 
valid for double pancake coils. 

2.1.5 Discretisation of the superconducting domain 

Due to their high aspect ratio, ReBCO superconductors are often approximated 
as 1-D objects, where the variation of the electromagnetic quantities along the 
thickness is neglected. However, in certain cases, this approximation is not 
correct. In particular, if we consider the cases shown in this paper, this happens 
for two reasons. First, for the uncoupled case (or a single tape) and low angles 
and low fields, the AC loss is dominated by the current penetration across the 
thickness of the tape (figures 151 and [TU|l . Second, for a field-dependent Jc and 
the uncoupled case (or a single tape), taking only one element in the thickness 
(or a ID approximation) neglects the influence of the local parallel field on 
Jc- This is because the parallel component of the field is anti-symmetric with 
respect to the mid-plane of the tape, so both its value at the mid-plane and 
its average in the thickness vanish. Taking only one element in the thickness 
induces significant errors for the uncoupled case and single tapes at applied 
fields below the self-field when Jc is not strongly anisotropic, as is the case of 
ReBCO coated conductors. 

The simulations use the following number of elements. For the coupled case, 
there are 1 element in the thickness and between 100 and 500 elements in the 
width per tape, with larger values for lower applied magnetic fields. For the 
uncoupled case and single tapes, there are between 1 and 20 elements in the 
thickness and 100 and 500 elements in the width. Higher numbers of elements 
in the thickness and the width are for lower applied magnetic fields. For the 
particular case of 0° (parallel applied field) and the uncoupled situation, we use 
50 elements in the thickness and 10 in the width. 
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2.2 Angular and field dependence 



In this article, we study the effect of three different field and angular depen- 
dences of the critical current density, Jc{B,9): isotropic superconductor, ellipti- 
cal anisotropy and a realistic field and angle dependence. This latter dependence 
presents a peak in both the ab and c directions and contains three different con- 
tributions with elliptical anisotropy (figure HJi). In the following sections, we 
study the error in the computations as a consequence of choosing a simplistic 
angular dependence, i.e. isotropic superconductor or elliptical anisotropy. In 
order to do so, we assume that the superconductor is perfectly described by the 
realistic description. Then, we simulate a measurement of the in-field critical 
current for an applied field in the ab and c directions (figure [^la). This is the 
data available in many experiments. With this data, we extract the parameters 
assuming an elliptical anisotropy. This simplification describes well the critical 
current in the ab and c directions but not for intermediate angles at large ap- 
plied magnetic fields. If we simplify further, to an isotropic superconductor, we 
could only expect to fit the critical current at one orientation. Moreover, it is 
not possible to fit the field dependence at low applied fields and high applied 
fields at the same time because the self-field changes the orientation of the local 
magnetic field. In this article we choose to take the parameters that fit the 
critical current at high applied fields in the ab direction. 



2.2.1 Realistic field and angle dependence 

Many ReBCO coated conductor tapes present an angular dependence with a 
peak in the ab and c directions. Reference |19| found that an expression based 
on three contributions that describes a coated conductor manufactured by Su- 
perpower, Inc. 113. Here, we use a simplified version of that field and angle 
dependence of Jc, which reproduces the main features of the experimental crit- 
ical current in |19) . 

JciB,0) = max{J,,ab[Bfab{e)lJcABfc{0)lJcABnm} (1) 

with 



(1 + BjBoab)^ ' 



JcabiB) = , (2) 



and 



fabiO) = \ cos^e + ul^sin'e, (5) 



fc{e) = Ju^cos^e + sin^e, (6) 



Me) = ■y/cos2e' + wf sin^e*. (7) 

We take the following values for the parameters: Jop = 4.9 • 10^° A/m^, Joi = 
3.2-101" A/m^, Boab = 4.6 mT, Boc 2.0 mT, Bq, = 32 niT, (3 ^ 0.48, a = 0.9, 
Uab = 8.3, Uc = 1.8 and Ui — 1.7. 
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2.2.2 Elliptical anisotropy 



Now, we consider the Blatter's [IH] scaling law for the field and angular depen- 
dence, also known as elliptical anisotropy: 

j,{B,e) = j,,,\BU{e)] (8) 

with 

"^^'^^^^ " (f + b/Bq,y ' 



/e(6') = ^COa'^e + ulsiT^e. (10) 

For this field and angle dependence we set the parameters in the following way. 
We take the parameters for which the critical current Ic as a function of the 
applied magnetic field Ba for the elliptical dependence, equations ([5|)- (fTU)) . fits 
the best to Ic for the realistic dependence, equations (H])-®, if the applied field 
is in the ab or c directions (see figure [2}3)- For low applied magnetic fields, 
the critical current is smaller than the integration of Jc{Ba,0) over the volume 
because the self-field increases the local magnetic field. The parameters for the 
elliptical anisotropy are Joe = 4.602 • 10^°A/m^, i?oe = 4.6mT, = 2.015 and 
/? takes the same value as for equations ([2]) and ©, /3 = 0.48. 



2.2.3 Isotropic superconductor 

Finally, the simplest case is that of an isotropic superconductor, with no angular 
dependence of Jc- 

'^^^^o{B,e) = j^-^^^. (11) 

For this field dependence, we take the same parameters as the field dependence 
for the elliptic case in the ab direction, equation ([3]): i?o = 4.6mT, Jq = 4.9 • 
10^" A/m^ and /3 = 0.48 (if we wanted to take the dependence in the c direction, 
the isotropic field dependence would be with Bq = Boe/ue = 2.283mT). 

With the isotropic dependence we do not find the parameters that fit the 
critical current for the realistic dependence. This is because an isotropic critical 
current density cannot describe the in-field critical current of a tape with an 
anisotropic superconductor, not even for a particular orientation. 



3 Results and discussion 
3.1 AC loss for a single tape 

In this section we present and discuss the AC loss in a single tape for an applied 
magnetic field with several orientations and the three angular dependences of 
section (results in figures Actually, we present the loss factor relative to 
the applied magnetic field amplitude, Q /B^ {Q is the loss per cycle and length 
and Bm is the applied field amplitude) or to the component of the applied 
magnetic field perpendicular to the tape surface, Q/B^^^. 

The results from both simulation methods agree (figures [3] and S]). However, 
there is a small discrepancy at large amplitudes because of the finite flux-creep 
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exponent for the FEM simulations (exponent 35). Indeed, when we increase the 
flux creep exponent to 101, this discrepancy vanishes (figure [3]). 

As expected, the AC loss decreases with decreasing the orientation angle 
[501 [TUl \T2\ [551 [23 ■ This is because the AC loss in a thin tape is dominated by 
the perpendicular component of the applied magnetic field and this component 
of the applied magnetic field decreases with decreasing the angle. Moreover, 
with decreasing the angle, the peak in the loss factor moves to higher applied 
fields because the perpendicular component of the applied field decreases. 

The shoulder at low applied magnetic fields for low angles (curve for 3° in 
figures [3| and [4]) is because of the current penetration across the thickness [14| . 
When the perpendicular component of the applied field is low, the AC loss due 
to the parallel component becomes important. 

The loss factor relative to the perpendicular applied magnetic field, Q/B^^^, 
is generally not independent on the angle (figure [1]), in contrast to the con- 
clusions of earlier work (TUl [H [221 [HI [301 [211 [22 • However, there actually is 
this angle independence for large and moderate angles (>30°). Apart from the 
shoulder at low amplitudes, for low angles (<15°) the peak of Q/B^^j. increases, 
shifts to lower fields and becomes slightly narrower. If we compare 3° with 90°, 
the peak for 3° is at around 3 times lower Bp^, resulting in around 3 times 
larger ac loss for i?por close to the peak for 3°. The difference in AC loss is simi- 
lar for large Bp^r- This angle dependence of Q/Bp^j. is caused by the anisotropy 
and field dependence of the critical current density. 

The results for the isotropic field dependence and the elliptical one reveal 
the different effects of the anisotropy and the field dependence, as follows. 

The increase of the peak in Q/Bp^^, its shift to lower angles and its narrowing 
when decreasing the angle are due to the field dependence because these effects 
are already present for the isotropic Jc (figure [It). In particular, the shift 
of the peak of Q/B^^^ to lower Bper with decreasing the angle is because of 
the following. For the same 5per, lower angles result in higher applied field 
magnitudes and, consequently, in a lower average Jc in the AC cycle. Then, the 
tape saturates for lower Bpor and the loss factor shifts accordingly. The increase 
in the peak of Q/B^^^ and its narrowing is due to the stronger field dependence 
of Jc as a function of Bper for lower Bp^r [5T]. Indeed, the field dependence of 
PT|) as a function of Bpcr is Jc ~ '^7(1 + -Bpor/(So| sin6'|))'^, resulting in an 
effective field constant i?o,off = Bq \ sin 6*1 which decreases with the angle 9 and, 
therefore, strengthens the field dependence with decreasing the angle. 

An elliptical anisotropy weakens the angular dependence of Q/B^^^. with 
respect to the reahstic or isotropic dependence (figure [U^). This is because 
the angular dependence in Jc of ([TO)) becomes fe{9) ^ Uc \ sin 6\ for large enough 
angles 9, resulting in a Uper dependence independent on the angle, Jc = Joe/(l + 
BperUe/ Boe)^ ■ The angle range where Jc becomes independent of i?pcr becomes 
wider for larger anisotropics, that is larger Ug. This explains the published angle 
independence of Q/Bp^^. as a function of Bp^ for Bi2223, with a large anisotropy 

m- 

The AC loss for the realistic case presents the main features of the isotropic 
case regarding the position, height and width of the peak. This is because the 
anisotropy of Jc is not very strong. The main difference from the isotropic case is 
that the angular dependence of Q/^por loss is smaller, especially for Bper below 
the peak. This is because for low fields, the realistic Jc approaches an elliptical 
dependence (figure [2^), therefore a weaker angular dependence of Q/B^^^.. 
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3.2 Magnetic field and current distribution in a Roebel 
cable 

The magnetic field distribution for the coupled and uncoupled cases are very 
different (figure For the coupled case, the superconductor shields the applied 
magnetic field as much as possible in all the cable volume. For the uncoupled 
case, the superconductor only shields the volume of each tape individually. This 
is in accordance to the published results for perpendicular applied magnetic 
fields [231 [3T]. In addition, for the uncoupled case, the superconductor shields 
the perpendicular component of the applied magnetic field within each stack of 
tapes but not in the gap between the stacks, where it concentrates. As a result, 
the angle of the magnetic field in this gap increases. Another consequence is that 
the magnetic field within each stack is parallel to the tape surfaces. Moreover, 
the magnetic field in between the tapes of the stack is uniform. Actually, this 
is always the case when the magnetic field in the horizontal separation between 
tapes is parallel to their surface f21| . 

In some cases it is not enough to take only one element in the thickness of the 
tape (ID approximation), see section [2.1.51 This is the case for the uncoupled 
situation and low applied magnetic fields. An example is for an applied magnetic 
field at 7° and 20 niT of amplitude, figure [51 where the AC loss is dominated 
by the current penetration across the thickness of the tape (figure [TUj) . On 
the contrary, for an applied field of, for example, 15° and 50 mT of amplitude 
(figure [6]) , one element in the thickness is sufficient for AC loss and magnetic 
field calculations. This is because the main AC loss contribution is from the 
penetration across the width of the tape and the average parallel magnetic field 
in the tapes thickness is nonzero. 

As expected; the current distribution between the coupled and uncoupled 
cases is very different (figures [6] and [8]). For the coupled case, the current 
penetrates roughly as in a monoblock. With an oblique applied magnetic field, 
the current penetrates faster from the top-left and bottom-right corners, figures 
[5^ and [5^. This is consistent with the calculations for one single tape [3] and the 
field penetration in figure [S^. For the uncoupled case, the net current in each 
tape is zero with qualitatively similar current penetration in each tape, except in 
the tapes at the boundaries of the cable (figures[6}3 and[8jD). The only similarity 
between the coupled and uncoupled cases is that the current distribution is anti- 
symmetric with respect to the central point of the tape. This is because of the 
geometry of the cable and the applied magnetic field. 

The current distribution for the ID approximation allows an accurate quan- 
titative comparison between the simulation methods (MMEV and FEM), rather 
than a qualitative comparison between current distributions from colormaps like 
figures [HI and [51 or the comparison of the AC loss in a logarithmic scale, such 
as in figures [5I4I9I and [TUl The simulation methods agree very well, better than 
for a constant Jc [H]. Moreover, the effect of the smooth E{J) relation for the 
FEM simulations is not very important for low applied field amplitudes [2Tj . 
With increasing the amplitude, the local electro-motive force due to the applied 
magnetic field increases, and so does the local electric field and J for a smooth 
E{J) relation. This is indirectly seen in the higher AC loss at high amplitudes 
for the FEM model (figures [51 and [TOl). 

The sheet current density, K, at the peak of the AC cycle presents the 
following features (figure [7|. Thanks to the inversion point-symmetry, it is 
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enough to study one half of the cable, for example the leftmost half in figure 
[51 For both cases, the sharp peaks in figure [7] correspond to the boundary 
between the regions with sheet current density equal and below its critical value 
{Kc ~ Jed). This is because at this boundary \K\ = Kc and the magnetic 
field is minimum (it vanishes for the coupled case, while for the uncoupled case 
only its perpendicular component vanishes). At the region in between the sharp 
peaks and the edge of the tape, \K\ becomes Kc- There, \K\ decreases toward 
the edge because the magnetic field increases. For the FEM simulations, the 
kinks and small peaks close to the sharp peaks are not necessarily numerical 
errors. A superconductor with a smooth E{J) relation presents a retarded 
response comparing to the critical-state model [52]. Then, at the peak of the 
AC field, the maximum current penetration is still not fully developed. For the 
uncoupled case, the sheet current density around the mid-width of the tapes 
presents a plateau (figure [TJj), as for perpendicular applied fields, transport 
currents and pancake coils [HJ [31]. The plateau is due to the difference in 
the local magnetic field at both surfaces of the tape and its height does not 
depend on Jc, as detailed in [3T]. This difference in field appears because of 
the inhomogeneous field created by the other tapes. Another issue is that the 
sheet current density close to the lowest y-edge of each tape is higher than for 
the highest y-edge. The reason is that at these places the sheet critical current 
density takes the critical value and the magnetic field concentrates at the top-left 
and bottom right corners of each stack of tapes, decreasing Jc- For the coupled 
case, the maximum sheet current density is higher because the superconductor 
shields better the magnetic field, decreasing the field and increasing Jc- The 
sheet current density in each layer of tapes is qualitatively different. In the first 
layer (black line in figure [Tji), the sheet current density roughly decreases with 
increasing y because the magnetic field increases (figure [5]). In the second layer 
(blue line in figure [7^), there is a double peak in the tape at higher y (a double 
peak for MMEV and single peak for FEM). This is because the magnetic field 
vanishes in the small region in between the peaks. 

When the simulations take into account the penetration into the thickness 
of the tape, they show that for the uncoupled case and low angles there is an 
significant penetration in the thickness (figure [S}d). Actually, for the situation in 
figure [Hb, the tapes are completely penetrated with current of both signs, with 
only a small difference corresponding to the height of the plateau in figure [7}d. 
For the coupled case, the distance of penetration across the thickness is much 
smaller than across the width, therefore the description with a sheet current 
density is correct in terms of AC loss (note that the actual aspect ratio of the 
tape is 4000). Moreover, the parallel magnetic field keeps the same sign in the 
thickness of the tape and, therefore, its average in the thickness is representative 
for the average critical current density. 

3.3 AC loss for a Roebel cable 

Again, the AC loss calculated from both simulation methods agrees with each 
other (figures [Ql fTT]) . There is a small discrepancy at high amplitudes or low 
angles and low amplitudes caused by the different physical model of the super- 
conductor (sharp and smooth E{J) relation for MMEV and FEM, respectively). 

The results for the coupled case describe the situation of untransposcd tapes. 
Then, the difference between the AC loss for the coupled and uncoupled cases 
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reveals the maximum possible reduction of AC loss by the transposition in a 
Roebel cable. For large applied magnetic fields or low angles, the AC loss for 
the uncoupled case is lower than for the coupled one (figures and ITU)) ) . The 
decrease by uncoupling the tapes is especially important for low angles, with a 
decrease up to three orders of magnitude for parallel fields and high amplitudes. 
For angles below 15° this reduction is still substantial for all the amplitudes. 
At high applied fields and angles of 7° or lower the loss reduces more than one 
order of magnitude compared to the coupled one. 

The AC loss for the coupled case presents the following features (figure ^ . 
It decreases with decreasing the angle until it saturates to the value for 0°. 
This decrease with decreasing the angle is less pronounced than for a single 
tape (figure [3]). The reason is that the coupled case behaves as a single block 
[131 [231 m] and the aspect ratio of the Roebel cable is relatively small, around 5. 
Then, the difference in the penetration length across the thickness and the width 
of the cable is less important than for the tape, resulting in a less pronounced 
angular dependence. 

For the uncoupled case (figure [TU]), the AC loss is qualitatively similar to a 
single tape (figure [3]). However, there are the following differences. First, the 
peak in the loss factor is at higher amplitudes because of the stacking effect 
[251 [31] [23] ■ Second, the peak (or peak at higher amplitudes) is wider than 
for a single tape because of the concentration of the perpendicular field in the 
gap between the stacks (figure [5]) [231 [H] . Finally, the peak at low applied 
fields and angles, due to the parallel applied field, becomes more visible. This is 
because two reasons. First, the peak due to the perpendicular applied field shifts 
to higher amplitudes, so both peaks do not overlap. Second, the contribution 
from the parallel field becomes larger than that of the perpendicular one. The 
cause is that the stacking effect reduces the AC loss per volume created by the 
perpendicular field but the loss created by the parallel field remains the same. 
For 0°, we can compare the simulations with the slab Bean model for a constant 
Jc [S31 113] taking Jc as the self-field critical current of the tape divided by its 
volume. The simulations are consistent with those Bean-model calculations 
(figure [TO)l . The higher and narrower peak of the loss factor for the simulations 
are due to the field dependence of Jc |SI]. Then, the Bean model is useful to 
determine whether the AC loss due to the parallel field is important, assessing 
the validity of the ID approximation at low angles without the need of lengthy 
simulations with several elements in the thickness. 

The loss factor relative to the perpendicular applied field Q/B^^^ (figure [TT]), 
does not depend only on Bper, as is the case for a single tape (figure [3]). With 
decreasing the angle, the shift in the peak and its increase (or the peak at high 
amplitudes for low angles) is less pronounced than for the single tape. However, 
for low i?por, the increase of the loss factor due to the parallel component of the 
applied field is more important. 

In real Roebel cables, the angular dependence of Q/B^^^ could be larger. 
This is because there can be misalignment between the tapes, they can be 
slightly deformed (like making a circular arch) or the superconductor layer could 
have a non- negligible roughness - see for example pictures in [27]. All these 
effects will result in larger contributions of the parallel applied field. In addition, 
coupling currents will also increase the AC loss from the parallel applied field. 

The AC loss for low angles will be important for long solenoids, like in 
some transformers. Therefore, we recommend to characterise the AC loss in 
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Roebel cables and tapes also at low angles of the applied field, therefore avoiding 
the error committed extrapolating Q/B^^^ for low angles assuming an angular 
independence of Q/Bp^^. 

4 Summary and conclusions 

In this article, we have presented the main features of the magnetic field, current 
distribution and AC loss for an applied magnetic field with arbitrary orientation, 
based on two independent numerical simulations. For the Roebel cable, we have 
taken into account the coupled and uncoupled cases, corresponding to the limits 
of very high coupling currents and negligible coupling currents, respectively. 
The simulations for the uncoupled case are also valid for a double pancake coil 
with no transport current, therefore the results and conclusions for the Roebel 
cable in the uncoupled case are also applicable to double pancake coils. The 
simulations have taken into account a realistic anisotropic field dependence of 
Jc- This dependence is in accordance with Jc measurements within around 20% 
error. Therefore, the results from the simulations are representative for real 
tapes, in contrast to those from simple dependences, such as the elliptical or 
isotropic ones. 

About the two simulation methods in this article, we have seen that all 
the results agree with each other, serving as a mutual check. In addition, the 
qualitative results for the field and current distributions are well explained by 
magnetostatic considerations. 

For both single tapes and Roebel cables in the coupled case, we found that 
the AC loss does not only depend on the perpendicular component of the applied 
field, opposite to the published works [IHl [HI [28l [29l [3ll |32]. For low 
angles between the applied field and the tapes plane and moderate and large 
amplitudes, the extrapolated AC loss from purely perpendicular applied fields 
assuming independence of the parallel component of the applied field is around 
3 times lower than the actual one, therefore being this difference not negligible. 
The discrepancy is even larger for low amplitudes. The AC loss also depends on 
the parallel component of the applied field for two reasons. First, at low applied 
fields and low angles, the AC loss due to the penetration across the thickness 
is important. Second, because for YBCO the field anisotropy in Jc is relatively 
weak. Then, Jc reduces significantly with increasing the parallel component of 
the applied field for the same perpendicular component of the applied field. 

For coupled Roebel cables, the AC loss due to the parallel and perpendicular 
components of the applied field is of the same order of magnitude for all the 
amplitudes. This is because of the low aspect ratio of the cable, around 5. 

Finally, we have found that the highest potential of Roebel cables is to 
reduce the AC loss at applied fields with low angles with the tapes surface. 
For angles of 7° or lower, the Robel cable reduces the AC loss more than one 
order of magnitude, comparing to a cable made of untransposed tapes with the 
same dimensions. The reduction can be of up to three orders of magnitude for 
perfectly parallel applied fields. This is a much better reduction than the factor 
2 for perpendicular applied fields [201 [H] ■ 

In conclusion, transposed cables, such as Roebel cables, should be used in 
windings with an important parallel applied field, such as the low- voltage wind- 
ing of transformers. For the characterisation of the AC loss of these cables, it 



12 



is necessary to measure also at low angles with the tape surface (< 15°). This 
consideration is also valid for double (or single) pancake coils to be part of a 
long solenoid, such as the high-voltage winding of a transformer. As future 
work, we propose comparison with experiments, the prediction of the influence 
of intermediate coupling currents and the role of striated tapes in the Roebel 
cable. 
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Figure 1: The sketch shows the origin of the angle 6 between magnetic field B 
and the surface of the tape. 
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Figure 2: In the top graph, the apphcd magnetic field is, from top to bottom: 
20, 60, 100 and 200 mT. In the bottom graph, the lines are for 0° and 90°, from 
top to bottom. 
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Figure 3: AC loss per unit length and cycle, Q, for a single tape with the 
realistic anisotropic field dependence of Jc in section 12.21 The angle of the 
applied magnetic field is 90, 30, 15, 7, 3 degrees in the arrow direction. 
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Figure 4: Loss factor Q/B^^^ as a fui3A;ion of the perpendicular applied field 
amplitude -Bpcr for a single tape for the anisotropic field dependences in section 
2.21 realistic, elliptic and isotropic (from top to bottom). The symbols distin- 
guish different angles of the applied field and the solid/dash lines distinguish 
the results from MMEV and FEM. 



Figure 5: The magnetic flux lines for the coupled case (up) are very different 
from the uncoupled one (bottom). These flux lines are calculated with the 
MMEV model at the peak of the AC cycle for 6=15° and B,„=50 mT. 
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Figure 6: This cross-section shows that the current distribution for the coupled 
(up) and uncoupled (bottom) cases are evidently different (the current distribu- 
tion is calculated with the MMEV model) . The situation is for the peak of the 
AC cycle for an applied field of 61=15° and 50 niT of amplitude. This current 
distribution is for the ID approximation: only one element in the superconduc- 
tor thickness. For better visualisation, the horizontal and vertical axis are not 
on scale and the represented tape thickness in not in scale with the horizontal 
separation between tapes. 
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Figure 7: The average current density J in the thickness of the tapes for both 
simulation methods agree (sohd and dash hues are for the MMEV and FEM 
resuhs, respectively). The the different lines (black, blue, green and red) are 
from the left border to the central layers of tapes in figures [5] and jS) The 
situation is the same as in those figures. 
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Figure 8: This cross-section shows that the current distribution for the coupled 
(up) and uncoupled (bottom) cases are evidently different (the current distri- 
bution is calculated with the MMEV model). The situation is for the peak of 
the AC cycle for an applied field of 9^7° and 20 mT of amplitude. This current 
distribution calculated with 20 elements in the thickness shows the details of 
the current penetration across the tapes thickness. For better visualisation, the 
horizontal and vertical axis are not on scale and the represented tape thickness 
is not in scale with the horizontal separation between tapes. 
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Figure 9: Loss factor, B„i is the applied field amplitude, for a Roebel cable 
in the coupled case. Solid lines with symbols are calculated by MMEV and 
dash lines with symbols are calculated by FEM. The angle of the applied field 
amplitude is 90, 60, 30, 15, 7, 3 and degrees in the arrow direction. 
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Figure 10: Loss factor, Bm is the applied field amplitude, for a Roebel cable 
in the uncoupled case. Solid lines with symbols are calculated by MMEV and 
dash lines with symbols are calculated by FEM. The angle of the applied field 
amplitude is 90, 60, 30, 15, 7, 3 and degrees from top to bottom. The black 
solid line is for the Bean slab model with a parallel applied magnetic field [S3] 
(see text). 
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Figure 11: The normalised loss factor relative to the perpendicular applied field 
amplitude, Q/Bpcx, for the uncoupled case does not only depend on -Bpcr only. 
Solid lines with symbols are calculated by MMEV and dash lines with symbols 
are calculated by FEM. The angle of the applied field amplitude is 90, 60, 30, 
15, 7 and 3 degrees from top to bottom in the arrow direction. 
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